Stability analysis of an autonomous system in loop quantum cosmology 
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We discuss the stability properties of an autonomous system in loop quantum cosmology. The 
system is described by a self-interacting scalar field (j) with positive potential V, coupled with a 
barotropic fluid in the Universe. With F = VV" /V''^ considered as a function of A = V' /V, 
the autonomous system is extended from three dimensions to four dimensions. We find that the 
dynamic behaviors of a subset, not all, of the fixed points are independent of the form of the 
potential. Considering the higher-order derivatives of the potential, we get an infinite-dimensional 
autonomous system which can describe the dynamical behavior of the scalar field with more general 
potential. We find that there is just one scalar-field-dominated scaling solution in the loop quantum 
cosmology scenario. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The scalar field plays an important role in modern cos- 
mology. Indeed, scalar-field cosmological models are of 
great importance in the study of the early Universe, es- 
pecially in the investigation of inflation. The dynam- 
ical properties of a scalar fields also make an interest- 
ing research topic for modern cosmological studies [ll, Q ■ 
The dynamical behavior of scalar field coupled with a 
barotropic fluid in a spatially flat Friedmann-Robertson- 
Walker universe has been studied by many authors (see 
P,[3,ll, and the first section of [1). 

The phase-plane analysis of the cosmological au- 
tonomous system is an useful method for studying the 
dynamical behavior of a scalar field. One always consid- 
ers the dynamical behavior of a scalar field with an expo- 
nential potential in the classical cosmology 0-01 or modi- 
fied cosmology [1,0. And, if one considers the dynamical 
behavior of a scalar field coupled with a barotropic fluid, 
the exponential potential is also the first choice [lol - [l3| . 
The exponential potential V leads to the fact that the 
variables F = VV" /V^ equal 1 and that A = V /V is 
also a constant. Then the autonomous system is always 
two dimensional in classical cosmology [5|, and three di- 
mensional in loop quantum cosmology (LQC) Q. Al- 
though one can also consider a more complex case with 
A being a dynamically changing quantity [l], [13, [l^ , the 
fixed point is not a real one, and this method is not exact. 
Recently, Zhou et al [H, [l^ introduced a new method 
by which one can make F a general function of A. Then 
the autonomous system is extended from two dimensions 
to three dimensions in classical cosmology. They found 
that this method can help investigate many quintessence 
models with different potentials. The goal of this paper 
is to extend this method for studying the dynamical be- 
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havior of a scalar field with a general potential coupled 
with a barotropic fiuid in LQC. 

LQC [ll,|il is a canonical quantization of homoge- 
neous spacetime based on the techniques used in loop 
quantum gravity (LQC) [13, [2l|. Owing to the homo- 
geneity and isotropy of the spacetime, the phase space of 
LQC is simpler than that of LQG. For example, the con- 
nection is determined by a single parameter c and the 
triad is determined by p. Recently, it has been shown 
that the loop quantum effects can be very well described 
by an effective modified Friedmann dynamics. Two cor- 
rections of the effective LQC are always considered: the 
inverse volume correction and the holonomy correction. 
These modifications lead to many interestin g res ults: the 
big bang can be replaced by the big bounce [22l | , the sin- 
gularity can be avoided j23l|. the infiation can be more 
likely to occur (e.g., see [24n28| ). and more. But the in- 
verse volume modification suffers from gauge dependence 
which cannot be cured and thus yields unphysical effects. 
In the effective LQC based on the holonomy modification, 
the Friedmann equation adds a — term, in which 
K = SttG, to the right-hand side of the standard Fried- 
mann equation 29]. Since this correction comes with a 
negative sign, the Hubble parameter H , and then d will 
vanish when p = pc, and the quantum bounce occurs. 
Moreover, for a universe with a large scalar factor, the 
inverse volume modification to the Friedmann equation 
can be neglected and only the holonomy modification is 
important. 

Based on the holonomy modification, the dynamical 
behavior of dark energy has recently been investigated 
by many authors f^, [sO, [3l| . The attractor behavior of 
the scalar field in LQC has also been studied [1^ Ull . It 
was found that the dynamical properties of dark-energy 
models in LQC are significantly different from those in 
classical cosmology. In this paper, we examine the back- 
ground dynamics of LQC dominated by a scalar field with 
a general positive potential coupled with a barotropic 
fluid. By considering F as a function of A, we inves- 
tigate scalar flelds with different potentials. Since the 
Friedmann equation has been modified by the quantum 
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effect, the dynamical system will be very different from 
the one in classical cosmology, e.g., the number of dimen- 
sions of autonomous system will change to four in LQC. 
It must be pointed out that this method cannot be used 
to describe the dynamical behavior of scalar field with ar- 
bitrary potential. To overcome this problem, therefore, 
we should consider an infinite-dimensional autonomous 
system. 

The paper is organized as follows. In Sec. [TTl we 
present the basic equations and the four dimensional dy- 
namical system, and in Sec. IIIIl we discuss the properties 
of this system. In Sec. IIV[ we discuss the autonomous 
system in greater detail, as well as an infinite-dimensional 
autonomous system. We conclude the paper in the last 
section. The Appendix contains the analysis of the dy- 
namical properties of one of the fixed points, P3. 



II. BASIC EQUATIONS 

We focus on the flat Friedmann- Robertson- Walker cos- 
mology. The modified Friedmann equation in the effec- 
tive LQC with holonomy correction can be written as 



1 



p_ 

Pc 



(1) 



in which p is the total energy density and the natural 
unit K = 8ttG = 1 is adopted for simplicity. We consider 
a self-interacting scalar field (f) with a positive potential 
V{(j)) coupled with a barotropic fluid. Then the total 
energy density can be written as p ^ p,p + p^, with the 
energy density of scalar field — ^(jP + V{(I>) and the en- 
ergy density of barotropic fluid p^ . We consider that the 
energy momenta of this fleld to be covariant conserved. 
Then one has 



(t) + 3H<j) + V' = 0, 
p^ + S'-fHp^ = 0, 



(2) 
(3) 



where 7 is an adiabatic index and satisfles = (7 — 
with being the pressure of the barotropic fluid, and the 
prime denotesthe differentiation with respect to the fleld 
(j). Differentiating Eq. ([1]) and using Eqs. 1^ and ([3]), 
one can obtain 



H 



1 



2(P7 + 



Pc 



(4) 



Equations (H])-® and ©-(HI) characterize a closed sys- 
tem which can determine the cosmic behavior. To an- 
alyze the dynamical behavior of the Universe, one can 
further introduce the following variables [1, Q : 



V 



Pc 



v 



where the z term is a special variable in LQC [see Eq. 
([1])]. In the LQC scenario, the total energy density p 



should be less than or equal to the critical energy density 
Pc, and thus < z < 1. Notice that, in the classical 
region, z — for p <^ pc- Using these new variables, one 
can obtain 



Pi 
H 



1 



1 - z 
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(6) 
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Using the new variables ([5]), and considering Eqs. ^ 
and (O, one can rewrite Eqs. (HJ-© in the following 
forms: 



dx 
dN 

dy 
dN 

dz 

In 

dX 

dN 



^,x-fxy^ + x 
x(l-2z). 



(8) 



^/6 



Xxy + y 



x(l-2z), (9) 
-37Z - 3z (1 - z) {2x^ - jx^ - 7y2) , (10) 

V6X^x{T~l), (11) 



where = In a and 



r ; 



VV" 
yi2 



(12) 



Note that the potential V{<j)) is positive in this paper, but 
one can also discuss a negative potential. Just as 33] has 
shown, the negative scalar potential could slow down the 
growth of the scale factor and cause the Universe to be in 
a collapsing phase. The dynamical behavior of the scalar 
fleld with the positive and negative potential in brane 
cosmology has been discussed by [3|. In this paper we 
are concerned only with an expanding universe, and both 
the Hubble parameter and the potential are positive. 

Differentiating A with respect to the scalar field 0, we 
obtain the relationship between A and F, 



dx- 



(13) 



If we only consider a special case of the potential, like 
exponential potential |5l4l3{. then A and F are both con- 
stants. In this case, the four dimensional dynamical sys- 
tem, Eqs. ([5|)- pT|) . reduces to a 3-dimensional one, since 
A is a constant. (In the classical dynamical system, the 
z term does not exist, and then the dynamical system is 
reduced from three dimensions to two dimensions.) The 
cost of this simplification is that the potential of the field 
is restricted. Recently, Zhou et al (Tfl [i3 considered the 
potential parameter F as a function of another poten- 
tial parameter A, which enables one to study the fixed 
points for a large number of potentials. We will follow 
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this method in this section and the sections that follow 
to discuss the dynamical behavior of the scalar field in 
the LQC scenario, and we have 



r(A) = /(A) + i. 



(14) 



In this case, Eq. p4|) can cover many scalar potentials. 

For completeness' sake, we briefly review the discussion 
of [13] in the following. From Eq. ([T3|) . one can obtain 



dX 



dV_ 
'V' 



(15) 



Integrating out A = A(y), one has the following differen- 
tial equation of potential 



dV 
V\{V) 



(16) 



Then, Eqs. and ((TE)) provide a route for obtain- 

ing the potential V = V{4)). If we consider a concrete 
form of the potential (e.g., an exponential potential), the 
dynamical system is specialized (e.g., the dynamical sys- 
tem is reduced to three dimensions if one considers the 
exponential potential for dX/dN = 0). These special- 
ized dynamical systems are too special if one hopes to 
distinguish the fixed points that are the common prop- 
erties of scalar field from those that are just related to 
the special potentials [l3| ■ If we consider a more general 
A, then we can get the more general stability properties 
of scalar field in the LQC scenario. We will continue the 



discussion of this topic in Sec. IIVI In this case, Eq. pT|) 
becomes 



(17) 



Hereafter, Eqs. ©-(HOl) along with Eq. (HZl) definitely 
describe a dynamical system. We will discuss the stabil- 
ity of this system in the following section. 



III. PROPERTIES OF THE AUTONOMOUS 
SYSTEM 



Obviously, the terms on the right-hand side of Eqs. 
(|5))- ([TU1) and pT)) only depend on x, y, z. A, but not on N 
or other variables. Such a dynamical system is usually 
called an autonomous system. For simplicity, we define 
= Fiix,y,z,X) = F,,^ ^ F2ix,y,z,X) = F2,^ ^ 
F3{x,y,z,X) = F3, and |^ = Fi{x,y,z,X) = F4. The 



dN 

fixed points {xcVct Zc, Xc) satisfy F^ = 0,i = 1,2,3,4. 
From Eq. p7|) . it is straightforward to see that x — 
0, A = or /(A) = can make Fi{x,y,z,X) — 0. Also, 
we must consider A^/(A) — 0. Just as [17] argued, it 
is possible that X^f{X) ^ and |^ 7^ when A = 0. 
Thus the necessary conditionfor the existence of the fixed 
points with a; ^ is A^/(A) = 0. Taking into account 
these factors, we can easily obtain all the fixed points of 
the autonomous system described by Eqs. (|8|l-(fT0|) and 
P?)) . and they are shown in Tab. I. 



TABLE I: The stability analysis of an autonomous system in LQC. The system is described by a self-interacting scalar field 
(j) with positive potential V coupled with a barotropic fluid p-y. Explanation of the symbols used in this table: Pi denotes the 
fixed points located in the four dimensionsal phase space, which are earmarked by the coordinates {xc,yc,Zc,Xc). A* means 



that A can be any value. \a is the value that makes /(A) = 0. M 



points. /i(A) 



dfW 



with A = 0, Xa. A 



2/(A„) + \a 



dX 



means the inverted matrix of the eigenvalues of the fixed 
U stands for unstable, and S stands for stable. 



Fixed-points 


Xc 


yc 


Zc 




Eigenvalues 


Stability 


Pi 














M^ = (0,-37,17,-3+17) 


U, for all 7 


P2 











A. 


M^ = (0,17,-37,-3+17) 


U, for all 7 














S, for7 = l,/i(0)>0 


P3 





1 








M^ = (-3,-37,0,0) 


U, for7= |, if /i(0) /O 
S, for7=|, if/i(0) = 


P4 


1 











M^ = (0,-6,0,6-37) 


U,for all 7 


Ps 


-1 











M^ = (0,-6,0,6-37) 


U,for all 7 


Pe 











A, 


M^ = (0,17,-37,-3+17) 


U,for all 7 


P7 


1 








A„ M^ 


= (-6, 6 - 37, l^/6Aa + 3, ^/QXaA) 


U,for all 7 


Ps 


-1 








A, M^ = 


(-6, 6 - 37, '\V&\a + 3, -^/GXaA) 


U,for all 7 


P.) 









Aa M^ = ( 


-Aa, —3 + |Aa, \i — 37, — Aa — fl{\a)) 


S, for /i(Aa) > Xa, and \a < 
U, for fi(Xa) < Aa and/or A^ > 3 


Pio 


fa 1 

V 2 A„ 1 







A„ 


See the Eq. (fT9)l 


U, for all 7 
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The properties of each fixed points are determined by 
the eigenvalues of the Jacobi matrix 



M = 



dx 


dy 


dz 


d\ 


dFj. 


dF2_ 


dF2 


dF2 


dx 


dy 


dz 


OX 


dFk 


OFj 


dFa 


dFs 


dx 


dy 


dz 


OX 


OFj 


dFj 


dFi 


dFi 


dx 


dy 


dz 


OX 



(18) 



{Xa,ya,Za,Xa) 



According to Lyapunov's linearization method, the sta- 
bility of a linearized system is determined b y th e eigen- 
values of the matrix Ai (see Chapter 3 of [341 )• If all 
of the eigenvalues are strictly in the left-half complex 
plane, then the autonomous system is stable. If at least 
one eigenvalue is strictly in the right-half complex plane, 
then the system is unstable. If all of the eigenvalues are 
in the left-half complex plane, but at least one of them 
is on the iui axis, then one cannot conclude anything 
definite about the stability from the linear approxima- 
tion. By examining the eigenvalues of the matrix M for 
each fixed point shown in Table I, we find that points 
Pi, 2, 4-8, 10 are unstable and point Pg is stable only un- 
der some conditions. We cannot determine the stability 
properties of P3 from the eigenvalues, and we will give 
the full analysis of P3 in the Appendix. 
Some remarks on Tab. I: 

1. Apparently, points P2 and Pg have the same eigen- 
values, and the difference between these two points 
is just on the value of A. As noted in the caption 
of Table I, A* means that A can be any value, and 
Aa is just the value that makes /(A) = 0. Obvi- 
ously, Xa is just a special value of A,, and point 
Pe is a special case of point P2. Pq is connected 
with /(A), but P2 is not. From now on, we do not 
consider separately the special case of Pg when we 
discuss the property of P2. Hence the value of Aa 
is contained in our discussion of A*. 

2. It is straightforward to check that, if Xc = Xc — 0, 
Uc can be any value ?/* when it is greater than or 
equal 1. But, if > 1, then Zc ^ 1 — < 
1, and this means that the fixed point is located 
in the quantum-dominated regions. Although the 
stability of this point in the quantum regions may 
depend on /(A), it is not necessary to analyze its 
dynamical properties, since it does not have any 
physical meanings. The reason is the following: If 
the Universe is stable, it will not evolve to today's 
pictures. If the Universe is unstable, it will always 
be unstable. We will just focus on point P3 staying 
in the classical regions. Then t/c = = 1, Zc = 
1 — l/y* = 0, i.e., for the classical cosmology region, 

P < Pc- 

3. Since the adiabatic index 7 satisfies < 7 < 2 (in 
particular, for radiation 7 = | and for dust 7 = 1), 



all the terms that contain 7 should not change sign. 
A more general situation of 7 is < 7 < 2 [351 . 
If 7 = or 7 — 2, the eigenvalues correspond- 
ing to points Pi, 2, 4, 5, 9 will have some zero elements 
and some negative ones. To analyze the stability 
of these points, we need to resort to other more 
complex methods, just as we do in the Appendix 
for the dynamical properties of point P3. In this 
paper, we just consider the barotropic fluid which 
includes radiation and dust, and 7 7^ 0, 2. Notice 
that if one considers 7 = 0, the barotropic fluid 
describes the dark energy. This is an interesting 
topic, but will not be considered here for the sake 
of simplicity. 

4. -V6 < Aa < \/6, Aa 7^ should hold for point Pg, 
hence -3 + ^Xl < 0. 

5. Aa > should hold, since j/c > for point Piq. The 
eigenvalue of this point is 



M = 



-37 
-3AQ7/i(Aa) 



V(2-7)(A?(2-7) + 87 + 2472) 



V-i + |7 - if: V(2 - 7)(A2(2 - 7) + 87 + 2472)7 



(19) 



Since we just consider < 7 
per, it is easy to check that (2 
87 + 247^) > is always satisfied. 

point is unstable with /i(Aa) 



< 2 in this pa- 

- 7)(A^(2 - 7) + 
And this 

be- 



dfjX) 
dX 



X=X 
3,3 
2 4 



+ 57 + 



ing either negative or positive, since 
_4r V(2-7)(A^(2-7) + 87 + 2472) is always pos- 
itive. 



Based on Table I and the related remarks above, we 
have the folloing conclusions: 

1. Points Pi, 2: The related critical values, eigenvalues 
and stability properties do not depend on the spe- 
cific form of the potential, since Ac = or A can be 
any value A*. 



2. Point P3: The related stability properties depend 
x=o' 



on /i(0) 



3. Points P4,5: The related eigenvalues and stability 
properties do not depend on the form of the poten- 
tial, but the critical values of these points should 
satisfy X^f{X) — since Xc ^ 0. 



4. Point Pg: It is a special case of P2, but /(Aa) 
should be satisfied. 







5. Points Pt %: Same as Pg, they would not exist if 

/(Aa) ^ 0. 
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6. Point P9,io: /(A^) = should hold. The fixed val- 
ues and the eigenvalues of these two points depend 

on /i(Aa) = ^ 



X = Xa 

Thus, only points Pi 2 are independent of /(A). 

Comparing the fixed points in LQC and the ones in 
classical cosmology (see the Table I of [13), we can 
see that, even though the values of the coordinates 
(a^c. He, ^c) are the same, the stability properties are very 
different. This is reasonable, because the quantum mod- 
ification is considered, and the autonomous system in 
the LQC scenario is very different from the one in the 
classical scenario, e.g., the autonomous system is four di- 
mensional in LQC but three dimensional in the classical 
scenario. Notice that all of the fixed points lie in the 
classical regions, and therefore the coordinates of fixed 
points remain the same from classical to LQC, which we 
also pointed out in an earlier paper [sTj . 

Now we focus on the late time attractors: point P3 
under the conditions of 7 = l,/i(0) > and 7 = 
4/3, /i(0) = 0, and point Pg under the conditions of 
A^ < 6,/i(Aa) > Aa,Aa < 37. Obviously, these points 
are scalar-field dominated, since — i/^(l/(l — Zc) — 
~ Vc) = 0- Foi" point P3, the effective adiabatic in- 
dex 70 — (p0 -\-P4>)/ P4> = 0, which means that the scalar 
field is an effective cosmological constant. For point Pg, 
l<t> — ^al'^- This describes a scaling solution that, as the 
universe evolves, the kinetic energy and the potential en- 
ergy of the scalar field scale together. And we can see 
that there is not any barotropic fiuid coupled with the 
scalar-field-dominated scaling solution. This is different 
from the dynamical behavior of scalar field with exponen- 
tial potential V = Vbexp(— Ak0) in classical cosmology 
[a-tlSf. and also is different from the properties of the 
scalar field in brane cosmology [s^l, in which A = const, 
(notice that the definition of A in 3] is different from the 
one in this paper) and F is a function of L{p{a)) and 
\V\. In these models, the Universe may enter a stage 
dominated by a scalar field coupled with fiuid when A, 7 
satisfy some conditions [1, HI . 

We discuss the dynamical behavior of the scalar field 
by considering F as a function of A in this and the pre- 
ceding sections. But F can not always be treated as a 
function of A. We need to consider a more general au- 
tonomous system, which we will introduce in the next 
section. 



IV. FURTHER DISCUSSION OF THE 
AUTONOMOUS SYSTEM 



The dynamical behavior of the scalar field has been 
discussed by many authors (e.g., see [ll. Isl-fisj). If one 
wants to get the potentials that yield the cosmological 
scaling solutions beyond the exponential potential, one 
can add a ^ term into the autonomous system (36| . 
All of these methods deal with special cases of the dy- 
namical behavior of scalar fields in backgrounds of some 



specific forms. By considering F as a function of A, one 
can treat potentials of more general forms and get the 
common fixed points of the general potential, as shown 
m and in the two preceding sections. However, 

as is discussed in [17], sometimes F is not a function of A, 
and then the dynamical behaviors of the scalar fields dis- 
cussed above are still not general in the strict sense. For 
a more general discussion, we must consider the higher- 
order derivatives of the potential. We define 



(i)r 



VV3 



(2)r 



17^' 



(3)r 



y/2 



(n)p _ 



VVr 



n+2 



in which Vn — 

dT 

dN 
d((i)F) _ 

dN 
d((^)F) _ 

dN 
d((3)F) 



y/2 ' 

= 3,4,5, 



(20) 



Then we can get 



dN 
d ((")F) 



dN 



rA + (^)F-2AF2 , 

"(i)rA + (2)r-2AF ((i)f)" 
"(2)FA + (3)r-2AF ((2)r)" 
WF-2AF ((3)r)" 

VGx ("'FA-H^"+^^r-2AF((")F) ,(25) 



(3)rA- 



(21) 
(22) 
(23) 
(24) 



To discuss the dynamical behavior of scalar field with 
more general potential, e.g., when neither A nor F is con- 
stant, we need to consider a dynamical system described 
by Eqs. ©-(dT]) coupled with Eqs. (IlI])-(Il5l). It is easy 
to see that this dynamical system is also an autonomous 
one. We can discuss the values of the fixed points of 
this autonomous system. Considering Eq. (Ill[) . we can 
see that the values of fixed points should satisfy Xc — 0, 
Ac = 0, or Fc = 1. Then, we can get the fixed points of 
this infinite-dimensional autonomous system. 

1. If Xc = 0, considering Eqs. ((8t- (fT0|) . one can get 

(j/c, Zc, Ac) = (0, 0, 0) or (j/c, Zc, Ac) = (0, 0, A*), and 
Fc, '•"^Fc can be any values. 



2. If Ac = 0, considering Eqs. ((Sll- dTO)) . one can see 
that the fixed points of (x,y,z) are {xc,yc,Zc) — 
(0,?/,,l-l/y2)and {xc,yc,Zc) = (±1,0,0). If Xc = 
0, Fc and ^"^Fc can be any values, and if Xc = ±1, 
(")Fc = 0. 

3. If Fc ^ 1, considering Eqs. ©-(Uni), one 

can get that the fixed points of (x,2/, z,A) are 
{xc,yc,Zc,Xc) = (0,0,0, A*) and {xc,yc, z^, Xc) = 
(±1,0,0, A*)- And (")Fc should satisfy ("^Fc = A'^ 
There are other fixed points, which will be dis- 
cussed below. 

Based on the above analysis and Table I, one can find 
that points Pi_io are just special cases of the fixed points 
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of an infinite-dimensional autonomous systems. Consid- 
ering the definition of F (see Eq. (IT^ ). the simplest 
potential is an exponential potential when Tc — 1- The 
properties of these fixed points have been discussed by 
many authors If Xc = and j/c = 0, this corre- 

sponds to a fluid-dominated universe, which we do not 
consider here. If Xc = ±1, Fc = and '■"-'Fc — 0, we 
do not need to consider the F and the ^"•'F terms. Then 
the stability properties of these points are the same as 
points P4_5 in Table I, and there are unstable points. 
The last case is {Xc,yc, Zc, Xc) = (0,t/*,l — l/j/J,0) and 
F, ("^F can be any value. To analyze the dynamical prop- 
erties of this autonomous system, we need to consider 
the '•"^Fc terms. We will get an infinite series. In order 
to solve this infinite series, we must truncate it by set- 
ting a sufficiently high-order (^^)f to be a constant, for 
a positive integer M, so that d((*^)F) /dN = 0. Thus 
we can get an (M -I- 4)-dimensional autonomous system. 
One example is the quadratic potential V = ^m'^<j>'^ with 
some positive constant m that gives a five dimensional 
autonomous system, and another example is the Poly- 
nomial (concave) potential V = M^~"'(/)" [s^ that gives 
an (n -|- 3)-dimensional autonomous system. Following 
the method we used in the two preceding sections, we 
can get the dynamical behavior of such finite-dimensional 
systems. 

In the remainder of this section, we discuss whether 
this autonomous system has a scaling solution. 

If Xc = 0, then F^ ^ 0, ("^F^ ^ 0, and the sta- 
bility of the fixed points may depend on the trunca- 
tion. As an example, if we choose '^'F = 0, then we 
can get a six dimensional autonomous system. The 
eigenvalues for the fixed point (xc, j/c, -Zc, Ac, Fc, '■^■'Fc) = 
(0,0,0, Af„F,,(i)FO, where = or = A*, is 

M^ = (0,0, 0,^7, -37, -3 +^7). 

Obviously, this is an unstable point, and it has no 
scaling solution. The eigenvalues for the fixed point 
{xc, Vc, zc, Ac, Fc, (i)Fc) = (0, 1, 0, 0, F„ (i)FO is 

= (0, 0, 0, 0, -37, -3 - 37). 

According to the center manifold theorem (see Chapter 8 
of [13, or [111), there are two nonzero eigenvalues, and we 
need to reduce the dynamical system to two dimensions 
to get the stability properties of the autonomous system. 
This point may have scaling solution, but we need more 
complex mathematical method. We discuss this problem 
in another paper (40| . 

We discuss the last case. If Fc = 1, we can consider 
an exponential potential. Then the autonomous system 
is reduced to three dimensions. It is easy to check that 
the values {xec,yec, Zec) of the fixed points are just the 
values (xc, Uc Zc) of points Pg-io in Table I. We focus on 
the two special fixed points: 

Fi : (xccVccZcc) = (-A/V6, ^1- A76,0), 

F2 : {xccVec, Zee) = (-v^7/A, V37(2-7)/(2A2),0). 



Using Lyapunov's linearization method, we can find that 
F2 is unstable and Fi is stable if A < 37. It is easy 
to check that pj — H^[l/{1 — Zec) — xl^ — — 
when [xecTVecZec) = (-A/a/6, a/1 - A^/e, 0). From the 
above analysis, we find that there is just the scalar-field- 
dominated scaling solution when we consider the au- 
tonomous system to be described by a self-interacting 
scalar field coupled with a barotropic fiuid in the LQC 
scenario. 



V. CONCLUSIONS 

The aim of this paper is two-fold. We discuss the dy- 
namical behavior of scalar field in the LQC scenario fol- 
lowing [igI . [17} . To further analyze the dynamical prop- 
erties of scalar field with more general potential, we in- 
troduce an infinite-dimensional autonomous system. 

To discuss the dynamical properties of scalar field in 
the LQC scenario, we take F as a function of A, and 
extend the autonomous system from three dimensions 
to four dimensions. We find this extended autonomous 
system has more fixed points than the three dimensional 
one does. And we find that for some fixed points, the 
function /(A) affects either their values, e.g., for points 
P4-10, or their stability properties, e.g., for points Pa^g. 
In other words, the dynamical properties of these points 
depend on the specific form of the potential. But some 
other fixed points, e.g., points Pi^2,ai'e independent of 
the potential. The properties of these fixed points are 
satisfied by all scalar fields. We also find that there are 
two later time attractors, but the Universe is scalar-field 
dominated since = at these later time attractors. 

The method developed by [H, [l3| can describe the 
dynamical behavior of the scalar field with potential of 
a more general form than, for example, an exponential 
potential But it is not all-encompassing. If one 

wants to discuss the dynamical properties of a scalar field 
with an arbitrary potential, one needs to consider the 
higher-order derivatives of the potential V{4)). Hence the 
dynamical system will extend from four dimensions to 
infinite-dimensions. This infinite-dimensional dynamical 
system is still autonomic, but it is impossible to get all 
of its dynamical behavior unless one considers Fc = 1 
which just gives an exponential potential. If one wants 
to study as much as possible the dynamical properties 
of this infinite-dimensional autonomous system, one has 
to consider a truncation that sets (^^)F = Const., with 
Af above a certain positive integer. Then the infinite- 
dimensional system can be reduced to MA-A) dimensions. 
And we find that there is just the scalar-field-dominated 
scaling solution for this autonomous system. We only 
give out the basic properties of this infinite-dimensional 
autonomous system in this paper, and will continue the 
discussion in the paper in [4(| . 

We only get the scalar-field-dominated scaling solu- 
tions, whether we consider F as a function of A or con- 
sider the higher-order derivatives of the potential. This 
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conclusion is very different from the autonomous system 
which is just described by a scalar field with an expo- 
nential potential Q . This is an interesting problem that 
awaits further analysis. 



It is easy to get the eigenvalues and eigenvectors of 
Let A denote the matrix whose columns are the eigen- 
values, and S denote the matrix whose columns are the 
eigenvectors, and then we have 
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Appendix: The stability properties of the Point P3 

In Sec. mil we point out that it is impossible to get 
the stability properties of the fixed point if at least one 
of the eigenvalues of M is on the iuj axis with the rest 
being in the left-half complex plane. The fixed point P3 
is exactly such a point. In this appendix, we use the 
center manifold theorem (see Chapter 8 of [13] , or (ssj ) 
to get the condition for stability of P3. The coordinates of 
P3 are (0, 1, 0, 0) and the eigenvalues are (—3, —87, 0, 0). 
First, we transfer P3 to P3 {xc = 0, yc = y ^ 1 = 0, Zc = 
0, Ac = 0). In this case, Eqs. jUl-dTni) and ^ become 



-x'-iy + iyj^ (l-2z) 

^ = iV6Ax(y + l) + (y+l) 



+ -7 ((1 + z) 
(Al) 

3x^ + ^7 ((1 + z) 

(A2) 



dz 

— = -37z-3z(l-z)[2x2-7a;2-7(y+l)2],(A3) 
dX 



— = V6A2(/(0) + /i(0)A)x, 



(A4) 



where we have considered that the related variables 
(x,2/,z,A) are small around point {xc,yc, Zc, Xc) = 
(0, 0, 0, 0). Therefore the following Taylor series 





(-') 
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x/6 
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0\ 
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-37 







1 





1 
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1 




I ) 




Vo 





1 


0; 



(A6) 



With the help of S, we can transform Ai' into a block 
diagonal matrix 











o\ 





-37 




















vo 








0/ 




(A7) 



where all eigenvalues of Ai have negative real parts, and 
all eigenvalues of A2 have zero real parts. 
Now we change the variables to be 



(AS) 



Then, we can rewrite the autonomous system in the form 
of the new variables: 









fx + ^\\ 


Y 




y 




Z 


z 


A 


vv 






V - / 



/dX\ 

dN 


f-3 \ 




(X] 






dY 
dN 


-37 




Y 


+ 


G2 


dZ 
dN 


^ 




z 




G3 


\dN/ 


V 000) 











(A9) 

where Gi = Gi{X, Y, Z, A), {i = 1, 2, 3, 4) are functions of 
A, Y, Z, and A. It is easy to get d by substituting the 
transformations x ~ X — -^ Z, y = Y + z — X, X = Z 
into the R.H.S. of Eqs. ((Mt - ((M| . 

According to the center manifold theorem [slj, there 
exists a C°°-center manifold 



1 



= l + z 



1~ z 

/(A) = /(0) + /i(0)A-K--- , 

can be used, where /i(0) = ^^^^ I^^^q. 

We can get the Jacobi matrix Ai' of the dynamical 
system Eqs. (IXT|) - (IX4)) as 



M' = 



/-3 

-37 §7 



V y 



(A5) 



W,lc = {{X,Y,Z,X) : X = hiiZ,X),Y = h2iZ,X), 
/i,(0,0) = 0,A.(0,0) = 0} 

such that the dynamics of (jA9p can be restricted to the 
center manifold. Jhi is the Jacobi matrix of hi, and 
hi{Z,X),h2{Z,X) are 

hi{Z, A) = AiZ"^ + A2ZX + A^X^ + • • • , (AlO) 
h2{Z, A) = BiZ^ + B2ZX + SaP + • • • . (All) 

We just consider the quadratic forms of Z and A in this 
appendix. 
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Considering the center manifold theorem, we have 
dhi{Z,~X) dZ dhi{Z,X) dX 



dX 
dN 
dY 
dN 



dZ dN 
dh2{Z,X) dZ 
'dZ dN 



dX dN' 
dh2{Z,X) dX 
'dX m' 



(A12) 
(A13) 



Inserting the Eqs. (|JTO)) and ^KU) into dX/dN, dY/dN 
in Eq. (|X9t and Eqs. (|A12p - (|A13|) . and comparing the 
coefficients of dX/dN and dY/dN, we get 



^1 = 0, ^2-^, ^3=0, 



B2 = 0, B3 = -. 



^ 12' 



(A14) 



Then, the dynamics near the origin is governed by the 
following equations. 



dZ 
dN 
dX 
dN 



= -Z^X + -/Z^X 



■jX' 



(A15) 
(A16) 



We consider two different values of 7 to get the stability 
properties of this system. This is because a different 7 
will give a different dynamical systems. The first one to 
be considered is dust, which has 7=1. Then, we have 



dZ 
dN 
dX_ 

dN 



(A17) 
(A18) 



According to Lyapunov's theorem, we can define a Lya- 
punov function to analyze the stability properties of a 
dynamical system. Different dynamical systems have dif- 
ferent Lyapunov functions, and one dynamical system 
can also have different Lyapunov functions. But all the 
Lyapunov functions U should satisfy U{x.) > at the 
original point (Chapter 2 of ^3T]). Then we can define 



C/i = i(Z^+P). 
Using Eqs. (|A17|) and ()Al8p . wehave 



Iz^-x^ 



(A19) 



(A20) 



According to Lyapunovs stability theorems, the system 
is stable if /i(0) > 0. 

Now we turn to considering radiation, which has 7 = | . 
Eqs. (jAlsp and (|A16p become 



dZ 
dN 
dX 
dN 



-^Vi(O), 



-2A-' 



-Z^X. 



(A21) 
(A22) 



We need to consider three possible cases: (a) /i(0) 7^ 0, 
(b) /i(0) = 0,Z{N = 0) = 0, and (c) /i(0) = 0,Z{N = 
0) 7^ 0, since these three different cases will bring out 
three different dynamical systems. 

If /i(0) 7^ 0, the Lyapunov function can be defined as 



1 



1 + Z7(6A) + A2 



(A23) 



where A = /i(0) ff /(O) > 0, and A - -/i(0) ff /i(0) < 
0. Then one can get 



dU2 _ 12^2 [(Z2-P)2+5A4] 



dN 



[6A + 6AA2 + Z2]' 



> 0. 



(A24) 



Then this point is an unstable one. 

_If /i(0) = and Z{N = 0) = 0, Eq. (|M2t becomes 
dX/dN — ~'2X^. Defining Lyapunov function. 



then 



C/3 = 1 + A2, 



^ - -4A^ < 0- 
dN 



(A25) 



(A26) 



If /i(0) = and Z{N = 0) 7^ 0, one can get Z = C 
from Eq. (jA21[) . with a non-zero constant C. Equation 
(|A22|1 becomes 



The Lyapunov function can be defined as 



(A27) 



U4=[l--^X 



(A28) 



Then we have 
dU4 

In 



C4 



eP)' <o. 



(A29) 



Obviously, according to Lyapunovs stability theorem, 
this point is stable as long as /i(0) — 0, regardless of 
Z{N = 0) = or Z{N = 0) 7^ 0. 
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